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Abstract 

Recently, a quasi-orthogonal space-time block code (QSTBC) capable of achieving a significant 
fraction of the outage mutual information of a multiple-input-multiple output (MIMO) wireless com- 
munication system for the case of four transmit and one receive antennas was proposed. We generalize 
these results to ut ~ 2" transmit and an arbitrary number of receive antennas nj^. Furthermore, we 
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completely characterize the structure of the equivalent channel for the general case and show that for 
all riT = 2" and the eigenvectors of the equivalent channel are fixed and independent from the 
channel realization. Furthermore, the eigenvalues of the equivalent channel are independent identically 
distributed random variables each following a noncentral chi-square distribution with 4nfl degrees of 
freedom. Based on these important insights into the structure of the QSTBC, we derive an analytical 
lower bound for the fraction of outage probability achieved with QSTBC and show that this bound is 
tight for low signal-to-noise-ratios (SNR) values and also for increasing number of receive antennas. We 
also present an upper bound, which is tight for high SNR values and derive analytical expressions for 
the case of four transmit antennas. Finally, by utilizing the special structure of the QSTBC we propose 
a new transmit strategy, which decouples the signals transmitted from different antennas in order to 
detect the symbols separately with a linear ML-detector rather than joint detection, an up to now only 
known advantage of orthogonal space-time block codes (OSTBC). 

I. INTRODUCTION 

In recent years, the goal of providing high speed wireless data services has generated a great 
amount of interest among the research community. Recent information theoretic results have 
demonstrated that the ability of a system to support a high link quality and higher data rates in the 
presence of Rayleigh fading improves significantly with the use of multiple transmit and receive 
antennas [1], [2]. A very important aspect is always the availability of analytical expressions 
to describe the stochastic nature of the channel under consideration as given in [1], [3] for the 
MIMO channel. This offers an opportunity to obtain, e.g., closed-form analytical formulas for 
the ergodic capacity or the outage mutual information of such MIMO channels. E.g., in [4], the 
probability density function (pdf) of the random mutual information for independent identically 
distributed (i.i.d.) MIMO channels was derived in the form of the inverse Laplace transform and 
a Gaussian approximation of the pdf was presented. In [5], the impact of MIMO channel rank 
deficiency and spatial fading correlation on the mutual information was analyzed. Furthermore, 
for correlated channels, the optimal transmit strategy and the impact of correlation on the outage 
probability were derived in [6]. 

There has been considerable work on a variety of new codes and modulation signals, called 
space-time codes, in order to approach the huge capacity of such MIMO channels. The perfor- 
mance criteria of space-time codes were derived in [7], [8]. One scheme of particular interest 
is the Alamouti scheme [9] for two transmit antennas. Later on, [10] proposed more general 
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schemes referred to as orthogonal space-time block codes (OSTBC) with the same properties as 
the Alamouti scheme like, e.g., a remarkably simple maximum-likelihood decoding algorithm. 
Interestingly, the combination of OSTBC with a MIMO antenna system can be represented 
equivalently as a single-input-single-output (SISO) system, where the channel gain is equal to 
the Frobenius norm of the actual MIMO channel. The performance of orthogonal space-time 
block codes [11]-[13] with respect to mutual information was analyzed (among others) for the 
uncorrelated Rayleigh fading case in [14], [15] and for the more general case with different 
correlation scenarios and line of sight (LOS) components in [16]. OSTBC exploit multiple 
antennas at both the transmitter and receiver in order to obtain transmit and receive diversity 
and therefore increase the reliability of the system. With the knowledge of the stochastic nature 
of the resulting equivalent channel due to the employment of OSTBC in a MIMO system, the 
loss in mutual information of OSTBC in subject to transmission rate, number of receive antennas 
and channel rank was quantified in [14], whereas in [15] a comparison of OSTBC with a system 
applying beam-forming was presented. 

Unfortunately, the Alamouti space-time code for two transmit and one receive antennas is 
the only OSTBC, which, to the best of our knowledge, achieves the maximum possible mutual 
information of a MIMO system [1], since we can not construct an OSTBC with transmission 
rate equal one for more than two transmit antennas [10], [17]. Therefore, [18]-[20] designed a 
quasi-orthogonal space-time block code (QSTBC) with transmission rate one for four and eight 
transmit antennas. By properly choosing the signal constellations as done in [21]-[27], it is 
possible to improve the BER performance with ML-detection for the codes given in [18]-[20]. 
The BER performance of QSTBC with suboptimal detectors has been analyzed in [28], [29]. 

The performance of QSTBC with respect to outage mutual information (OMI) for the special 
case of one receive antenna and four or eight transmit antennas was analyzed via simulations 
in [19] and [30] and it was shown, that the QSTBC are capable to achieve a significant 
portion of the MIMO-OMI. Furthermore, it was shown in [19], that QSTBC in conjunction with 
optimal (nonlinear) and suboptimal (linear) detectors provide a tradeoff between performance 
and complexity. The key achievements of this paper are as follows 

• We generalize the results in [19] to 2" transmit and an arbitrary number of receive antennas. 

• We show, that due to the employment of QSTBC the eigenvalues of the resulting equivalent 
channel are pairwise independent and identical (i.i.d) noncentral chi-square distributed with 
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4?T,R degrees of freedom (xln^l'^nc)) with noncentrality parameter 5nc- Furthermore, we show 
that the eigenvectors of the equivalent channel are independent of each channel realization, 
i.e. they are constant. 

In other words, we first show that the combination of QSTBC with a MIMO system results also 
in a equivalent channel similar to OSTBC and then fully characterize the stochastic nature of 
this equivalent channel. Based on these important insights, we are able to provide the following 
results 

• an analytical lower bound for the outage probability achieved with QSTBC, which is tight for 
low signal-to-noise-ratios (SNR) values and also for increasing number of receive antennas; 

• an upper bound on the outage probability, which is tight for high SNR values. For the case of 
four transmit and an arbitrary number of receive antennas we derive analytical expressions 
for this bound. 

• finally, we exploit the special structure of the QSTBC and apply a new transmit strategy, 
which decouples the signals transmitted from different antennas in order to detect the 
symbols separately as in the case of OSTBC. The performance of this linear detector is 
equivalent to the non-linear maximum-likelihood (ML) -detector in [19]. 

The remainder of this paper is organized as follows. In Section |lll we introduce the system 
model and establish the notation. The design of QSTBC for 2" transmit antennas is shown 
in section [nil The complete characterization of the equivalent channel model and other key 
achievements of this paper are described in section |Wl As an possible application of the results 
in this section, we present the analysis of the outage probability achieved with QSTBC in 
section |Vl followed by some simulations and concluding remarks in section IV-BI and |Vll 

II. System model 

We consider a system with ut = 2" transmit and ur receive antennas. Our system model is 
defined by 

Y = G„^H + N, (1) 

where G„j, denotes the (T x ut) transmit matrix, Y = [yi,...,y„^] the (T x n^) receive 
matrix, H = [hi, ... , h.„^] the {ut x nj^) channel matrix, and N = [rii, . . . , n„^] the complex 
(T X nj^) white Gaussian noise (AWGN) matrix, respectively. An entry {nu} of N (1 < i < nji) 



February 1, 2008 



DRAFT 



5 



denotes the complex noise at the ith receiver for a given time instant t {1 < t < T). The real and 
imaginary parts of nu are independent and A/'(0,nr/(2SNR)) distributed. An entry of the channel 
matrix is represented by {hji} E hi and describes the complex gain of the channel between the 
jth transmit (1 < j < ut) and the ith receive (I < i < nj^) antenna, where the real and imaginary 
parts of the channel gains are independent and normal distributed random variables and hj is 
CJ\f{uii, I) distributed, where m, is the channel mean or Ricean component. The channel matrix 
is assumed to be constant for a block of T symbols and changes independently from block to 
block. The average power of the symbols transmitted from each antenna is normalized to one, 
so that the average power of the received signal at each receive antenna is ut and the signal- 
to-noise ratio (SNR) is p. It is further assumed that the transmitter has no CSI and the receiver 
has perfect CSI. 



III. Code construction 

A space-time block code is defined by its transmit matrix G„j,, which is a function of the 
vector X = [xi, . . . ,Xp]^. The rate i? of a space-time block code is defined as i? = p/T. In 
this paper, we focus on rate one QSTBC with length ut = T, therefore p = ut- Now, let us 
split the vector x into two vectors, Xodd and Xcvcn. for reasons that will be clear later on. The 
elements of x with odd index j are collected in Xodd and with even index in Xevcm respectively. 
Both parts of x are given as 



Xodd 



Si 



'"t/2 



rs~ x 



'"t/2+1 



(2) 



with Si, . . . , Snj, E C, where C C C denotes a complex modulation signal set with unit average 
power, e.g. M-PSK. Furthermore, T E C"^/^^"^/^ is a unitary matrix. More details on T and its 
effect on the detection scheme will be discussed in section IIV-GI 

Starting with the well known Alamouti scheme [9] for ut = 2 transmit antennas as a 



G2(a;i,a;2) 



Xi 



n 



X2 



-XT 



the generalization of the transmit matrix for the QSTBC with tit = 2" {tit > 4) is done in the 
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following recursive way 



G "T 

2 



G "T 

2 



Tlx 
T ■ 



- + 1 





where = xi,...,x„j, and the diagonal x "t/2 matrix ©„j, is given by 0„y 

diag 

Example 3.1: For the case of tit = 4 transmit antennas we have 



G,{{x,} 



j=i) 



Xi X2 X-s 



X4 
2^2 



2 ^2, 



In this work, we use the Alamouti scheme as the basis in order to construct the rate one QSTBC. 
However, it is also possible to construct QSTBC with rates lower than one based on other 
OSTBC [10], [12]. In the following section, we perform channel-matched filtering as the first 
stage of preprocessing at the receiver in order to obtain the equivalent channel model, followed 
by the decoupling of the system model in two parts. Afterwards, we analyze the eigenvalues 
and the eigenvectors of the resulting equivalent channel, leading to important insights of the 
properties of QSTBC. Noise pre-whitening as the second stage of preprocessing at the receiver 
is considered in section HV-FI 

IV. Signal Processing 

First of all, we briefly review the usual MIMO fading channel without any coordinated coding 
and the impact of OSTBC on the MIMO channel in order to provide a better insight into the 
properties of QSTBC. 



A. MIMO channel without any coordinated coding 

In this case, after channel matched filtering to ([l]), we have 

HH^ = VDDV^ , (3) 

where H = VDU^ is the singular value decomposition (SVD) of H, where the unitary matrices 
U,V contain the eigenvectors of H. The joint density function of the eigenvalues /ii, . . . , /i^ of 
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HH^ in DD in the Rayleigh fading case (rrij = 0) is given as [1], [3] 

1 



m\Kr, 



(4) 



where Km,n is a normalizing factor, n = max{nr, Jt-j?} and m = min{nT,nR}. It is obvious, 
that the eigenvalues are not independent of each other and it is well known that the matrix of 
eigenvectors V depend on the actual channel realization. 



B. Equivalent channel for OSTBC 

In case of OSTBC, the following holds for the transmit matrix 



Starting with ([T]), after some manipulations and channel matched filtering one arrives at 



where 



H 



nn riT 

EE 

i=i j=i 



E E \h,^\' 











riR riT 

EEI/^. 

i=ij=i 



(5) 



Since there is no interaction between the elements of x, the equation above can be decomposed 
into p parts. The resulting equivalent channel for each element of x of the OSTBC is then a 
single-input-single-output (SISO) channel given as 



(6) 



i=i j=i 

which is equal to the Frobenius norm of the actual MIMO channel matrix H. 

In case of the rate one QSTBC discussed in this paper, the actual MIMO channel is also 
transformed into a equivalent channel given as Hi^. Differently from the OSTBC the equivalent 
channel of QSTBC is still a MIMO channel, however with very interesting properties like 
constant eigenvectors and i.i.d. eigenvalues following a noncentral X4„^((5„c) -distribution as 
derived in the following. 
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C. Channel-Matched Filtering 

After rearranging and complex-conjugating some rows of Y the system equation in (HJ) can 
be rewritten as 

y' = H' x + n', (7) 



where H^^ = [(H^^ J^, . . . , (H^^^J^, . . . , (H;^_„J^]^ and H;^_, is given as 



(8) 

Thus it appears that the H^^ ^ ({^jj}j=i)' with [{hjiY-'L^ = hu, . . . , hnj,i, are obtained recur- 
sively, where the recursion starts with n-r = 2, 

H2,j = H2^j(/lij, /i2i) = 

— /)* /)* 

"^21 "'li 

In order to perform channel-matched filtering we multiply (H^^)^ from left to Q to get 

y =H„^x + n , 

where the noise vector n" = (H^^)^n' is spatially colored and H^^ is given as 

K^K + L^L 



(9) 



H 



K^K + L^L 



(10) 



where K = Hi^ ({/^ji} Jd) and L = Hi^ ({^ji}"^-r+i 



D. Decoupling of the system model 

An important property of the QSTBC the system in ^ can be decoupled into two parts due 
to the special structure of H"^ as described in the following. The decoupling comes from the 
fact that for QSTBC, it holds that [21] 



^nTi (-'^odd) ■ Gm^fXcvcn) ~l~ G„ (Xcvcn) ' G.fi„ (Xodd 

) = Vx 



TLX ^ 



(11) 



where Xodd = Xodd ® [ 1 = [xi, 0, X3, 0, . . . , 0]^ and Xcvcn = Xcvcn ® [ lY- The 

property in (fTTTl is very crucial, because this enables a simple maximum-likelihood decoding 
algorithm. Assuming perfect channel estimation is available, the receiver computes the following 
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decision metric over all possible transmit matrices and decides in favor of the transmit matrix 
that minimizes the following decision metric based on ([l]): 

|Y-G„,(x)-H||| = tr{(Y-G„,,(x)-H)^(Y-G„,(x)-H)} 

= tr{Y^Y- Y^G„,^(x)H 



(12) 



Gfnj, (Xj;v(;n)H -f- Ngven ; 



-(Y^G„,(x)H)^ +H^G„,(x)^G„,(x)H} . 

After some manipulations, we arrive at 

t^lY^j^Yodd + YQjjjG„j,(xodd)H + H G„^(xodd)Yodd + H G„^( ^odd) Gj^j, (Xodd )H+ 

^even^cvcn ~l~ Yg-^pjjGyjy (XpY^jj^H -|- H G^^^ (Xcvcn) Y^vgn -|- H G^^^ (XgYgjj) Gy^j, (x^YgQ^Hj- , 

where tr{-} is the trace function. Yodd and Ycvcn are given as 
Yodd = G„j,(xodd)H + Nodd and Y 

cvc 

respectively. The above decision metric can be decomposed into two parts, one of which 

ti'jY^dYodd + Y^<jG„j,(xodd)H + H^G^^(xodd)Yodd + H^G^^ (xodd)G„,j,(xodd)H} 
is only a function of G„j,(xodd)5 and the other one 

tr-{^ Ygygj^ Ygvcn Yg-^^gjjG^j, (xgygjj)H -(- H Gj^^(xgvgn)Ygvgjj -|- H G^^(xgvgjj)G^j,(xgvgjj)Hj' , 

is only a function of G„y(xevcn)- Thus the minimization of ([T2b is equivalent to minimizing 
these two parts separately. Note that, due to the processing at the receiver, the property in (fTTTl 
is projected on the channel matrix H^^ in ®. The decoupled parts depend either on Xodd or 
Xgven given in Q. 

Thus, it is now possible to write an decomposed system model for each part based on 
The decomposed system model for the part with Xodd (and similarly for Xgvcn) can be rewritten 
as 

Yodd = H^ixXodd + n . (13) 

For illustration, we present two examples for the case of n-r = 4 and ut = 8 transmit antennas. 
Example 4.1: (ut = 4 transmit antennas) In this case, H4 ^ in (EJ) is given as 

hii h^i h^i 

h* —h* h* 

'hi 'Mi "'Si 

hii —h2i 

h* 
"'li 



H' 



4,i 



—h* 
'hi 

-hsi 



-h* 



-h* 
'hi 



h* 
'hi 
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and H4 appears in ^ as 








iot2 












ai 





—ia2 




X2 


-ia2 





ai 












ia2 









X4 



n" , (14) 



where ai and 02 are given as 

«! = ^ ^ and 02 = ^ 2Im(/i*./i3i + hlJi2i), (15) 

j=i j=i 1=1 

respectively. From (fT4l) . it is now directly obvious, that the system equation can be decoupled 
into two parts, which then can be considered separately. For the case considered in this example, 
the decomposed system model for Xodd (and similarly for Xcvcn. cf. (fTSl) ) can be written as 



Yodd = H2 



Xi 



n 



(16) 



with a non-orthogonal 

tti ia2 
—ia2 Oil 

Example 4.2: (ut = 8 transmit antennas) The same procedure applied here results in a H4 
given as 

«! ia2 ict3 a4 
—ia^ —04 ai ia2 

^4 — ««3 — ««2 Oil 



= H4 

2 



where 



riR 8 



ai 



"3 



i=l j=l 1=1 

hl,ih2i + hl^Jiji + h*^^h4^j), and 



6i J ) 



(17) 



j=i 



04 = ^2Re(/iij/i7i + /i8j/i2i - hl^h^i - hlih^). 

i=l 

The general case of arbitrary ut = 2" and very important insights about the eigenvalue decom- 
position, the eigenvalues themselves and the eigenvectors of the equivalent channel Hiir , which 
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are crucial and necessary for further analysis, e.g., the derivations of the lower and upper bound, 
are provided in the following section. 



E. Eigenvalues and Eigenvectors of the equivalent channel model 

In order to completely characterize the equivalent channel with respect to the eigenvalues and 
eigenvectors, we first look at the properties of matrices with certain structures and then show 
that the equivalent channel matrices fulfill this special structure. Let the matrix Mat, where 
N 



^ = 2" ^, has the following recursive definition 



MAr(ai, . . . ,q;7v) = 

Similarly 

NAr(aAr+i, . . . , Q!2Ar) = 

where the recursion starts with 

M2(ai,a2) = 



Miv(ai, 



2 



.an) Njvfai 



2 

, Oat) 



Mivfai, 



> oin) 



, aN_] 



(18) 



Mjv (asiv , . . . , Q!2Af) Niv (ctAT+l, . . . , Q!3iv ) 



(19) 



—ia2 tti 



and N2(a3,a4) 



— ^4 ia^ 



Remark 4.1: The matrices M2 and N2 have the following eigenvalue decompositions 



M2 = V2S2Vf and N2 = V2T2Vf (20) 



where 



V2 



1 
71 



1 1 
—i i 



and 



S2 










ai + a2 








i4 







ai — a2 


' ul 





= i 


as — 04 















as + 04 



T2 {{ai}U) 

Immediately the following question follows: Is their any structure how to derive the eigenvalues 
of the matrices of higher A^, i.e., if the eigenvalues of Miv are given, how can we compute 
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the eigenvalues of Mjy. (Note that, if the eigenvalues of Mjv are given it is straightforward to 
derive the eigenvalues of Njv). In order to answer this question we are able to state the following 
lemma, where the arguments of Ma? and Nat are omitted. 

Lemma 4.1: Let Mat, Nat be as given in (IT8l),(IT9l), then Mat, with N = 2""\ra > 2 have 
the following eigenvalue decomposition: 

= V^S^V^ and Njv = Va^TjvV^ , (21) 

where 

= (h ® Yn] Un (In ® V2) , (22) 



iTN 

2 





l=N/2+l 



S JV 

2 



r 1^/2 





2 



Ml 



N/2+1 



(23) 



Tat — Tat {{(^i}'i=N+i) — 

({aj^fa^) + zS^ ({a^Sv+i) 

and 

[n^],^. = 5[2j-l-z]+(5 

with 5[-] denoting the delta function, giving 5\l\ = 1 for / = and 6[l] = for / ^ 0, and [riArjij 
denotes the -element of the N x N permutation matrix IlAr. 

Proof: The proof is given in Appendix U 
It is important to realize that Sn and Tat are constructed with different arguments. 

The from H"^ in (flOb in even and odd block structure re-sorted equivalent channel matrix 
H:^ in (fT3l) has exactly the same structure as Mat. Therefore, Lemma 14.11 can be directly 
applied to Hi^. To emphasis the usefuUness of the resulting property of the QSTBC, we are 
able to state the following theorem. 
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Theorem 4.1: The left and right eigenvectors of the equivalent channel in (flSt of QSTBC, 
which fulfill the recursive construction rule of (HJl are given by eq. (l22b and therefore constant 
for any arbitrary channel realization. 

Remark 4.2: Another important aspect of Lemma WA\ is the fact that the eigenvalues in S^v 
can be obtained simply by adding the eigenvalues of Siv and Tjv in an appropriate manner as 
done in (|42ll (cf. Appendix IJ), which will be used in the following analysis of the QSTBC. 

Lemma 4.2: Let 5*, u^^ be as in Lemma WJ\ The eigenvalues of the equivalent channel 
matrix of the QSTBC are given by the recursive equations (l23t . Let Siir = D^^Di^ and 

1 < j < the eigenvalues of S^r. Then for any ut and 



(f^D^ = ^^^ir' where fii 

riji, the eigenvalues {fli^^y of ^Di^rD:^ are obtained as follows 



2 2 

UR 



^hfAf;^h„ i<j< 



i=l 



(24) 



where the matrices A{ with j = 1, 3 . . . , — 1 and nx = 2",n > 2 are given as 



A] 



A^ 



A] 
-B' 



Bi 



~2' 



(25) 



with W„ 



2 • 



z0„^A^^ and j' 
Proof: The proof is given in Appendix HH 
Theorem 4.2 ( [31]): If hj is CA/'(mj,I) and P is an ut x ut matrix then hf^Ph/ has a 
noncentral xK^nc) distribution if and only if P is idempotent (P^ = P), in which case the 
degrees of freedom is A; = 2rk(P) = 2tr{P} (where rk(P) and tr{P} denote the rank and trace 
of P, respectively) and Snc = m^Prrij. 

Lemma 4.3: The matrices A;^^ are Hermitian and idempotent. 
Proof: The proof is given in Appendix ITTII 
From Lemma 14.21 and 14.31 it is now possible to prove the following theorem. 

Theorem 4.3: Let DuxDz^t defined as in Lemma l4~2l be the diagonal eigenvalue matrix of 

22 >^ >^ 

the equivalent channel matrix of an QSTBC, which fulfill the recursive construction rule of 
Then for any and nji, the eigenvalues of ^DiirDnr are pairwise independent and 



identical noncentral chi-square distributed with Aur degrees of freedom. 
Remark 4.3: It is important to realize that 0^^ = 0„y and 0^^A{^0„y 
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Remark 4.4: Recall from section iTVi that we have decoupled the system into two independent 
parts. The derivation above holds therefore for both parts, i.e., each eigenvalue appears twice, 
once for each part. 

Proof: The proof of Theorem 14.31 is given in Appendix |IV| 



F. Noise pre-whitening 

Since n in ([T3t is colored noise, the next step is to perform pre-whitening. With the knowledge 
of the theorems 14 . 1 1 and l4~3l it is easy to compute the pre-whitening filter Fpw at the receiver now. 
To this end, we need just the eigenvalue decomposition of Hiir given as = V 

2 2 2 2 2 

with SiiT = DiirDiiT. Therefore the pre-whitening filter is given as Fpw = D^V^. By 

2 2 2 2 2 

multiplying Fpw from the left to STJ[ we arrive at 



Yodd = Hxodd + w , 

where the entries of w are mutually i.i.d. Gaussian processes again. 

Example 4.3: In the case of = 4 transmit antennas H in (E^ is given as 



H 



Pi 



and 



-ifil 



and jll 



(26) 



(27) 



(28) 



G. Linear maximum likelihood detection 



From theorems 14.11 and 14.31 it is now possible to determine F adequately, resulting in a 
attractive system equation, which allows a very simple but effective ML-decoding. To emphasize 
this property we formulate the following corollary. 

Corollary 4.1: By choosing the matrix F in Q as F = Vi^r, (l26b can be rewritten as 



Yodd = E)-rS + W . 



(29) 



At this point, the elements of s~ (and also s"*") are completely decoupled, since they experience 
no interference from each other. Thus, a linear ML-detector is able to detect the symbols (or 
elements) transmitted from the antennas separately. 
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Proof: The matrix H in (l26b can be decomposed as 



H = D-tV 



H 



2 ~ 



(30) 



where D:^^ is a diagonal matrix containing the singular values of H. Since Vi^ is constant 
for all channel realizations, we can set T = Vi^ without any knowledge of the current channel 
realization. Using (l30b in (l26b results in (l29b . That concludes the proof. ■ 
Example 4.4: For nr = 4 transmit antennas, V2 and D2 are given as 

III 



1 1 







Example 4.5: For ut = 8 transmit antennas, we have the following V4 

1111 



-t —I 
i —i 
1 -1 



-i i 
1 -1 



and D4 = A/4diag(/ig, 



As 



. , jlj) with 

^(«l+02+«3-"4y4^ = -^(ai+a2-03+"4y, 
^(Ql-Q:2+a3+a4)/4^ aud /ig = V^(^ 



)/4. 



V. Performance analysis 



(31) 



Based on these new insights, we provide some performance analysis in this section, where 
we focus on the case of Rayleigh fading (nij = 0). 



A. Outage Probability Pout 

The mutual information of a MIMO system with ut transmit and nj^ receive antennas with 
no CSI at the transmitter and perfect CSI at the receiver by using the optimal transmit strategy 
is given as [1] 

/ = log2det fl„^ + ^HH 

'in this paper, we use the same terminology as in [1], i.e. we use the term capacity only in the Shannon sense and distinguish 
therefore between the concept of outage mutual information (OMI) and capacity. 
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The portion of the mutual information achieved with QSTBC is 

Iq = — log2 det f + Ad^^Di^) = — log2 n f 1 + ^-^{K.,A ■ (32) 

The outage probability Pout achievable with QSTBC is defined as the probability that Iq is 
smaller than a certain rate R, i.e. 

PoutiR,nT,nji,p) = Fi[Iq < R] . 

Unfortunately, the exact analysis of Pout is not available. Therefore, we provide a lower and 
upper bound in the following. 
1 ) Lower bound: 

Proposition 5.1: The outage probability Pout is lower bounded by 

Po„i(i?,nT,n^,p)>l-exp(-^(2^-l)j — k\ ^^^^ 

Proof: By using the arithmetic mean - geometric mean inequality, i.e. 



L , L 



Wa/"^ < y ^a;, a/ > 



L 

1=1 1=1 

with equality if and only if ai = 02 = ■ ■ ■ = we obtain an upper bound for Iq (and therefore 
a lower bound on Pout) given as 

log, (1 + -^ai J = , (34) 

where ai = X]r=i X]J=i l^i^l^ general case of arbitrary ut similar to the cases riy = 4 

and ut = 8. The lower bound on the outage probability Pout can be written as 




Pout{R,nT,nn,p) = Ft[Iq < R] > Pr[/^ < R] = Fi 



P 



Since ai is chi-square distributed random variable with 2nTnR degrees of freedom. Pout is given 
as [32, p. 310,3. 351(1)] in (l33t . That concludes the proof. ■ 
Corollary 5.1: The lower bound in (l33t gets tight for low SNR values or when ur increases. 
Proof: The proof is given in Appendix |Vl 
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2) Upper bound: Using the positive definiteness of Ti^T)^, (l32b can be lower bounded as 

iQ>—\ogAi+{— 

riT y \nT 
Thus, the upper bound on Pout is given as 

^ 




R 

For the special case of = 4 transmit antennas, Po^i is given as 

Since 4:{fi\y {jliy is a product of two chi-square distributed random variables, both with m = 
Aur degrees of freedom, Po^t is given by [33, p. 365, eq.9.9.34] 

■'• 1 Hf fln(|) + 2*(*+l) + sif5)(|)'" 

where ^ is the Psi function [32, p. 943, eq.8.360] and V is the Gamma function [32, p. XXXI]. 
Note that for high SNR a useful and simple approximation of the outage probability can be 
obtained by retaining only the first term (i.e. A; = 0) of the series within the upper bound. 

Some simulation results of the performance of QSTBC and their interpretation are presented 
in the following section. 

B. Simulations 

In Fig. m the OMI of QSTBC Iq with our new transmit strategy and our linear detector is 
compared with the nonlinear ML-detector and ZF-detector in [19]. Additionally, the OMI of a 
MIMO system with ut = A and ur = 1 h depicted. From the Fig., we observe that our new 
transmit strategy outperforms the ZF-detector of [19] and achieve the same portion of mutual 
information as the non-linear ML-detector presented in [19]. 

In Fig. 121 the performance of QSTBC in terms of OMI with ut = A and ut = 8 antennas is 
depicted for > 1. For ur = 1, the performance with = 8 is similar to the case of n-r = 4 
transmit antennas (depicted in Fig.EJ, i.e. we achieve a significant fraction of the OMI. However, 
by increasing the number of receive antennas, we observe in Fig. |2l that the performance of 
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6 




-10 -5 5 10 15 20 

SNR[dB] 



Fig. 1. 10% Outage mutual information (OMI) of a MIMO system, our new approach, and the ML-and ZF-detector from [19] 
with riT ~ 4: and ur — 1. 




Fig. 2. 10% Outage mutual information of a MIMO system and our new approach with ut ~ 4,nT = 8 transmit and ur > 1 
receive antennas. 
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QSTBC with tit = 8 as well as with = 4 is dramatically reduced in terms of achievable 
OMI. 

In Fig. 121 Pout of QSTBC with = 4 transmit and ur = 1 to Ufj = 3 and ur = 6 receive 
antennas is depicted. From the Fig. we observe that our lower bound on the performance of 
QSTBC with respect to Pout gets tight for increasing number of receive antennas. Even the 
upper bound performs very well and shows to be useful. The performance of QSTBC with 




SNR 

Fig. 3. Outage probabilities of QSTBC (dashed lines), upper bound (dotted lines), and lower bound(solid lines) for wt = 4 
transmit and different size of receive antennas ur, Rate=4. 

respect to Pout is depicted in Fig. |4| for ut = S transmit antennas. Similarly to the case of 
nr = 4 transmit antennas, the lower bound gets tight by using many antennas at the receiver 
side. 

VL Conclusion 

In this paper, we generalized the QSTBC in [19] to 2" transmit and an arbitrary number of 
receive antennas. Our most important results are given in the theorems 14 . 1 1 and |431 in section HV-EI 
and reveal the useful properties of the resulting equivalent channel for QSTBC. In more details. 
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SNR [dB] 

Fig. 4. Outage probabilities of QSTBC (dashed lines) and lower bound(solid lines) for ut = 8 transmit and different size of 
receive antennas nn, Rate=4. 

theorem 1431 state that the eigenvectors of the equivalent channel are independent of each channel 
realization, i.e. they are fixed. This property can be directly exploited for preprocessing, which 
allows a very efficient linear ML detection (Corollary 14.11) . In more details, we developed a new 
transmit strategy, which allows to apply a linear ML-detector as in the case of OSTBC, such 
that the symbols from different antennas can be detected separately. The performance of our 
linear detector in terms of mutual information is equal to the nonlinear ML-detector in [19]. 
The theorem 14.31 offers insights in the statistics of QSTBC. It is proved that the eigenvalues of 
the equivalent channel are i.i.d. noncentral xlnal^nc)- These insights were used to derive upper 
and lower bounds on the outage probability with QSTBC. It was shown, both analytically and 
via simulations, that the lower bound gets tight for increasing number of receive antennas and 
also in the low SNR-regime. From simulation results, we observed that the QSTBC approaches 
the outage mutual information only in the case of ri/j = 1 receive antenna. By increasing the 
number of receive antennas, the loss in terms of mutual information increases unbounded with 
the SNR. 
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Appendix I 
Proof of Lemma |4lT] 



In the following, the arguments of M^r, N^v, Sat and Tat are omitted occasionally in order to 
increase the readability of the paper. 

Proof: The proof is done by the principle of induction. We start with the initial case M4. 
From (flSl . it follows 



M4 



M2 N2 
-N2 M2 

V2 
V2 



V2S2Vf V2T2Vf 
-V2T2Vf V2S2Vf 



S2 T2 

^T2 S2 



Vf 

vf 



V2 

V2 

V2 

V2 



^2 ^2 











/"2 

1 1 
-^2 f^2 



n 



H 



Vf 

v,^ 



V2 
V2 



vf 
vf 



n 



H 



Vf 

vf 



V4 



with 114 given as 



1000 
0010 
0100 
0001 



S4 = diag(/i4,/i4,/i|,/i|), where 

nl = iul , ^1 = ^1 + iul ^ jj,l = iul , and //^ = /i^ + iul. 



which is equivalent to 



82(01, 02) — iT2(a3, 04) 

82(01, ^2) + ^T2(a3, 04) 



n 



H 



(35) 



(36) 



(37) 
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The same procedure applied to N4 results in a T4 given as T4 = diag(^'], z/|, . . . , uf) with 



• 4 J 4 4 I • 4 

?/i2 ) ^IIQ 2/4=1/2+ ^/^2- 



which is equivalent to 

T,{{ai}l,) = n4 



T2(«7, tts) - «S2(a5,«6) 

T2(a7,a8) +iS2(Q!5,a6) 

Now assume that the following hypothesis holds for = ^/2, i.e. 



n 



H 



Mk 

2 



V V? and N 



K_ 
2 



then the following inductive step is true 





Mk 

2 


Nk 

2 


IP 


2 


SaV? 


VaT 

2 


£ V K 

2 T 




-Nk 

2 


2 . 




-Va 

2 


TaV? 

2 — 


VaS 

2 


if V K 
2 — 




Yk 

2 







Sx 

2 


Ta 

2 




2 












2 . 




-Tx 

2 


Sa 

2 . 







* A 
2 - 





Va 

2 







Va 

2 



n 







2 







Q 



n 



v^ 











vf 

2 



(^l2 ® Va) Uk (Ik ® V2) Sk (Ik ® Vf 



v^ 



Va 



* A 



where 







fJ'K 




2 


2 






-Z^A 


/Ua 


2 


2 



2 



and = diag(/i]^, /i^, . . . , /if) with 

/-I - - / - 

= A ~ ^z^A and = ^]^+ iv 
22 2 

which is equivalent to 



n 



Sa(W,LJ-zTa(M 



V/ = 2,4,6, 







K 







S-(MzLi)+^TA({az};i, 



+1' 



(38) 



(39) 



(40) 



(41) 



(42) 



n 
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The same procedure applied to N/< results in the same V4 and = diag(i^^, uj^, . . . , i^^) 
with 

K+l K+l K+l K+l 

= - H^K and z/j^ = + i/i/ V/ = 2, 4, 6, . . . , (43) 

22 22 



which is again equivalent to 



Tk({«,}^5^^^^^) + zSk({«,}^+^,) 



E 



Since the initial case of = 4 is true and the inductive step is true, the statement in (OTTl is 
true for all N = 2"^. That concludes the proof. ■ 

Appendix II 
Proof of Lemma lOl 

Proof: The proof is done by the principle of induction. The outline of the proof is as 
follows. For the initial case of rir = 4, we need the eigenvalues {jliy and (^'2)^ = 2, i.e. 

the Alamouti scheme, as indicated in (l36b . The first step is therefore to construct the eigenvalues 
for riT = 4 with the eigenvalues of ri^ = 2. Using (l42b and (l43l) . we observe that the eigenvalues 
for n-r = 4 can be also obtained with the eigenvalues for ut = 8, which is the second step 
revealing an important instruction of constructing eigenvalues fj.K,^K from jJiK,^^. It follows 
the hypothesis and the inductive step concluding the proof. 

Now, we start with the well known Alamouti scheme. By applying the Alamouti scheme 
{ut = 2), — Di^DiiT as well as Xodd and Xcvcn are only scalars. Thus, the only eigenvalue of 

nx 2 2 

DiDi of the decomposed system model for the part with Xodd (and similar for Xcvcn) is given 
as 

riR njt nT=2 ur 

{il\f = J2Mhu,h,,) = E E l^^-^l' = E ^^^"-2. AX->2. (44) 

i=l i=l j=l 1=1 

where hfc^z,i = [hki, • • • , huf and 



= A2 = I2 (45) 



Similarly, 



'.~.1n2 



Mhh hii) = hf-.4,iA2h3^4,i (46) 



1=1 1=1 
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We are now able to start with the initial case of tit = 4. The first eigenvalue of the QSTBC 
for the part with Xodd (and similar for Xcvcn) is given as in (l42t (with (/i^j,)^ = and 

EH " "^ 1 

= {fi\{hu,...,hii)f = fLl-ivl = ^-{ai{hii, . . . ,h4^i) + a2{hii, . . . ,h4i)) (47) 



(El 



3i + ^Ii^2i)) 



j=l 



-/)* /)* /)* —h* 

"'3i "'ii 'hi "'2i 



■4,i. 



i=l 



■4,j< 



A2 

A2 

A2 

A2 



hi^4,i - i 



h* —h* —h* h* 
"-li 'hi "'3i 'Mi 



A2 
A2 



hl^4,i 

(48) 



•4,i, 



02 
-02 



A2 
A2 



hi- 



Z^"-l^A,ic) 



where 



A2 










«02A2 





A2 




-i02A2 







— 


1 


A2 - 


-i02A2 






2 


«02A2 


A2 



hl^4,i — hf^4jA4hi^4 



(49) 



In an analogous manner, we get (z^]) given as 
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On the other hand, with (l42b we have 

= if^dhii,...,hsi)) =2^^((/i8) +(/^8)) = 2^7; ^ 



i=l 



i=l 



j=l \j=l 



^li ^2i 









02 







-0. 



hi_.8,i 



Q A2 

A2 



A2 
A2 



»8,i 



Zl A^^-^^^i 



i=l 



A2 



A. 











02A2 
-02A2 







E>f- 



1=1 



riR . . 



A2 



A2 -i02A2 

202A2 A2 







02A2 
02A2 



hi- 











A2 -i02A2 

Z02A2 A2 



hi- 



»8,i 



i=l 








) 





^4 





hi^s,, = hf-.8,^^ (I2 ® Ay hi^s,, (50) 



i=l 



riR 



j=l 



A^ 







A\ 



h5^8,i 



nR 



2'^1^4,jA4hi^4^j + -h^g jA4h5^8,; 

(51) 



In an analogous manner, we get {ulY given as 



nR 



" 1 



i=l 







^0fcAi 



-iSkAl 



h5->8,i 
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Since the eigenvalues and (t'l)^ can be obtained very easily in a similar way, we omit the 
derivations here. 

Comparing (BTTl with (l49b shows that in order to get the eigenvalues of tt-t = 8, the eigenvalues 
of = 4 have to be expanded by using the Kronecker product of I2 and A4 , divided by | 
in order to incorporate the channel entries h^i, . . . , hsi as given in (l50b . Actually, this can also 
be observed in the expansion from tt-t = 2 to nr = 4 by comparing (l44b with the first addend 
in (EHll. 

Now assume that the following hypothesis holds 



i=l 



K 

Ai 



hi- 



>2kA 



Similarly 



i=l 



20feAi 







hi- 



»2fc,: 



then the following inductive step is true 



r,iJ+i-\2 
2fc J 



^(/i{(/iii, . . . , h2ki)f =F {^({.hu, . . . , /i2fci))' 



|53l,|5l 



i=l 



\^l^k,i^k+l-^2k,i. 



i=l 



"A-r 







hi^fc,j 









hfc-|-l_>2fc,i 



A;,i^^fc+l— ►2fc,; 





k^k 



^^[hiLfc.ihfcH 



k+1^2k,ii 







i^e.Af Af 



hi-^fc,i 

hfe-|-1^2fc,i 

hi^fc.i 

hfe-f-l-»2fc,i 



hf^2fc,i A2fc'^^hi^2fc,i - ^ h^ 



2fe,; 



j=l 



i=l 



TiQkM. 



±^QkAi Ai 



(52) 
(53) 



(54) 



(55) 



hl^2/c,i ! 

(56) 
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with j = 1, 3, . . . , — 1 and j' = That concludes the proof. 



Appendix III 
Proof of Lemma 1431 

Proof: With ESI (A^ )^ and {A{+^)^, j = 1, 3, . . . , "t/2 and f = are given as 



2 



(Any )^ 

2 



Thus, A;^';^"''^ are only Hermitian, if A-t^ is Hermitian. Since A2 is Hermitian, it follows that 

~2~ 

A:^^, nr = 2",j = 1,3,..., "t/2 are Hermitian. Similarly, 



4 



2A^ 

i;-£T. Jr\. ri'j' 



=Fi 2 A^^ A:^^ 

2 A-'' A-'' 

~2~ ~2~ 



(57) 



Thus, A{';^+^ are only idempotent, if A^^ is idempotent. Since A2 is idempotent, it follows that 
A^^, ut = 2",j = 1,3,..., "t/2 are idempotent. That concludes the proof. ■ 



Appendix IV 
Proof of theorem |43I 



Proof: We first prove the independency of the eigenvalues. Since the matrices An^ are 
Hermitian and idempotent, we are able to rewrite (l24l) as 



riB. 



■riR 



i=l 



i=l 



Independency between the eigenvalues in the Gaussian case is given if and only if the eigenvalues 
are uncorrelated, i.e. 



which is fulfilled if 



E[(A^„ h,)«A:; h,] = \fj,j^k 



(A^^J^AL = A^„ AL = Wj,jy^k. 



(58) 



By applying the eigenvalue decomposition to (l58t . one has to distinguish between the case, 
where the eigenvalues are given as 



(I - e„JAl^ (I + Bn^Ai^ , (I + Bn^Ai^ (1-0 



(59) 
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and 



fl 



fl 



21 ) 

2 



(60) 



with j' = /c' = and j' 7^ k'. Note that the entries on the /th-diagonals of the matrices 
A, where / = ±{1,3,. ..,72^ — 1} (/ = represents the main diagonal, / > above the main 
diagonal, and / < below the main diagonal), are equal to zero. Due to the special structure of 

the matrices A, it follows that (iTGnr)-^^!^ is orthogonal to (I±©„j,)Ai2: and (I=f©„j,) A-tj, 

2 22 

is orthogonal to (I =F 0nT)Any , with j' ^ k'. Thus, the eigenvalues in ( l59b and (l60b are zero 

2 

and therefore the eigenvalues in (l24l) are independent. 

The probability density function (pdf) of the eigenvalues can be obtained from (l24b 

as follows. The rank (rk(-)) of A2 is 2. Furthermore, 



rk(A/ 



2"J 



rk(UA^2"U) = rk 



[I + ©nj,) A2„_i 





rk((I 



rk((I 



© A-' 





© lA-^' 

= rk(A^l_i; 



where U contains the eigenvectors of A^n . Since rk(A2^ 
2, and thus the following holds 

^ h 


where V is a unitary matrix. With (l6n) . the pdfs are given as 



2, the matrices A:^ have all rank 



V^(A^,JV 



(61) 



nT\2\ 



p tr 



i=l 



p tr 



2=1 



H 



I2 




which is the sum of squares of 2nji independent complex normal distributed variables, i.e. a 
noncentral chi-square distribution with Aur degrees of freedom. That concludes the proof. ■ 



Appendix V 
Proof of corollary 15.11 



Proof: The inequality (0?!) is tight only, if the ratio of two eigenvalues, i.e. r 
for all i ^ j. From this it follows that it has to be shown that the following holds 

lim Pr(r = 1) = 1 . 



1, 
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Since the eigenvalues are chi-square distributed with each Anj^ degrees of freedom, the ratio of 
the eigenvalues is distributed as follows 

which is the well-known F distribution [33, p. 365, eq.9.9.35]. Therefore, when nji goes infinity, 
the F distribution is given as 

lim {h{r,nfi)) = 6{r — 1) , 

n^— ►oo 

where S is the delta distribution. It follows that the lower bound gets tight for increasing n^. The 
lower bound is also tight for low SNR values, which is obvious after rewriting (l32b as follows 

Furthermore, (1 + ^ai)^ ^ C for small SNR. As an example, ( = (^02)^ for ut = 4. 
Therefore, as the SNR gets smaller, the lower bound gets tighter. That concludes the proof. ■ 
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